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Abstract: 

We have calculated the complete 0{a) and 0{as) radiative corrections to the Higgs- 
production process e"'"e~ ttH in the Standard Model. This process is particularly 
interesting for the measurement of the top-quark Yukawa coupling at a future e"'"e~ col- 
lider. The calculation of the 0{a) corrections is described in some detail including, in 
particular, the treatment of the soft and collinear singularities. The discussion of numeri- 
cal results focuses on the total cross section as well as on angular and energy distributions 
of the outgoing particles. The electroweak corrections turn out to be sizable and can reach 
the order of ±10%. They result from cancellations between electromagnetic, fermionic, 
and weak bosonic corrections, each of which are of the order of ±10%. 
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1 Introduction 



One of the main tasks at future colliders is the unraveling of the mechanism of elec- 
troweak symmetry breaking. In the Standard Model (SM) electroweak symmetry breaking 
is implemented via the Higgs-Kibble mechanism, which is responsible for the generation 
of particle masses and which leads to the prediction of a physical scalar particle, the Higgs 
boson. The mass of the Higgs boson is expected to be in the range between the current 
lower experimental bound of 114.4 GeV [Q] and 1 TeV. Electroweak precision tests yield 
a 95% CL upper limit on the mass of the SM Higgs boson of 211 GeV A SM Higgs 
boson in the mass range up to 1 TeV will be discovered at the LHC, provided it exists and 
has no exotic properties. However, for the complete determination of its profile, including 
its couplings to fermions and gauge bosons, experiments in the clean environment of an 
e"'"e~ linear collider [IHj are indispensable. 

The fermion masses in the SM are generated from the Yukawa interactions via a finite 
vacuum expectation value v = {\/2G^)~^^'^ ~ 246 GeV of the Higgs-boson field. The 
strength of the coupling of the physical Higgs boson to fermions is fixed hj Qfj^ = ruf/v 
in lowest order. The measurement of these Yukawa couplings is therefore important in 
order to verify the mass-generation mechanism of the SM. By far the largest Yukawa 
coupling is the top-quark-Higgs-boson coupling {g^m ~ 0.5) owing to the large top-quark 
mass. Its measurement is therefore especially interesting. 

A direct access to the top-quark Yukawa coupling is provided by the process e+e^ 
ttH [ Ej if the Higgs-boson mass is not too large, i.e. Mh ~ 100-200 GeV. This process 
proceeds mainly through Higgs-boson emission off top quarks, while emission from inter- 
mediate Z bosons plays only a minor role, making it suitable for a determination of Quu- 
If the Higgs boson is light, i.e. Mh ~ 120 GeV, a precision of around 5% can be reached at 
an e+e~ linear collider operating at y/s = 800 GeV with a luminosity of / L dt ~ 1000 fb~^ 
[IHj. Combining the ttH channel with information from other Higgs-pro duct ion and decay 
processes an even better accuracy can be obtained in a combined fit [Ej- Furthermore, by 
investigating the process e"'"e~ — >• ttH bounds on non-standard physics in the top-quark 
Yukawa coupling can be derived [IHIII. 

Assuming an experimental precision of a few per cent, a thorough understanding of 
the background to the ttH final state is necessary [|H], and precise theoretical predictions 
for the signal process e"'"e~ — > ttH at the per-cent level are needed. This requires the in- 
clusion of radiative corrections, a rather complicated task for a process with three massive 
unstable particles in the final state. As a first step, the process e+e" — >• ttH can be treated 
in the approximation of stable top quarks and Higgs bosons. The corresponding leading- 
order total cross section is already known for a long time [E]. The 0{as) corrections to 
the total cross section within the SM have been calculated first in the "effective Higgs- 
boson approximation" valid only for small Higgs-boson masses and very high energies [ 
ITU] . A calculation of all diagrams for the dominant photon-exchange channel has been 
performed in Ref. [iHj, while the full set of diagrams has been evaluated in Ref. [H^. 
Within the Minimal Supersymmetric Standard Model (MSSM) the 0{as) corrections to 
the production cross sections of neutral Higgs bosons in association with heavy quark pairs 
(tt, bb) have been calculated in Ref. [E] for the photon-exchange channel. In Ref. [ ITIj 
all QCD diagrams have been taken into account, while the SUSY-QCD corrections have 
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been worked out in Ref. [E]- Furthermore, the 0{as) corrections to the Higgs-boson 
energy distribution both in the SM and the MSSM have been discussed in Ref. [il6j. 

Recently, considerable progress has been achieved in the calculation of the electroweak 
corrections to e"'"e~ — > ttH. Results for the electroweak 0{a) corrections in the SM have 
been presented in Refs. [ElE]- A calculation of some top- mass enhanced electroweak 
corrections in the SM, the MSSM, and the two-Higgs doublet model has been performed 
in Ref. 

We have accomplished a completely independent calculation of the 0{a) electroweak 
corrections to the process e"'"e~ ttH in the SM as well as of the 0{as) QCD corrections. 
First results for total cross sections have already been presented in Ref. [1201 including a 
comparison of our results with those of Refs. [ El ■ While we find good agreement with 
Ref. [CHI 5 results differ from those of Ref. [El at high centre-of-mass (CM) energies 
and close to threshold. Moreover, we reproduced the results of Ref. [[T2| for the QCD 
corrections to this process. 

In this paper we present details of our calculation. Single hard-photon radiation is 
included using the complete matrix element and combined with the virtual corrections ap- 
plying both phase-space-slicing and subtraction methods. Higher-order corrections from 
initial-state photon radiation (ISR) are taken into account in the leading-logarithmic ap- 
proximation. Moreover, we supplement our previous work with results for distributions 
in the energies and production angles of the Higgs boson and the top quark. Results on 
electroweak corrections to differential cross sections have not yet been presented in the 
literature before. 

The paper is organized as follows: in Section |2 the calculation of the virtual and real 
corrections is described. Section jS] contains a discussion of numerical results. The paper 
is summarized in Section EJ Further useful information on the calculation of the matrix 
elements is collected in the appendices. 

2 Calculation of radiative corrections 

2.1 Conventions and lowest-order cross section 

We consider the process 

e"(pi,cTi) +e+(p2,fT2) — ^ t(A;i,ri) +t(A;2,r2) + H(/i;3), (2.1) 

where the momenta pa, h of the particles and the helicities of the electrons, aa, and 
top quarks, Xj, are given in parentheses. The helicities take the values aa,Ti = ±1/2, 
but we often use only the sign to indicate the helicity. The electron mass is neglected 
whenever possible, i.e. it is kept finite only in the mass-singular logarithms related to 
ISR. This implies that the lowest-order and one- loop amplitudes vanish unless ai = — o"2. 
Therefore, we define a = ai = —a2- The particle momenta obey the mass-shell conditions 
Pi = P2 — ^1 — ^2 — ''^t 5 ^3 — ^H- For later use, the following set of kinematical 
invariants is defined: 

s = (p^+p^y^ 
Sij = {ki + kj) , i,j = 1, 2, 3, 

tai = iPa-k)^ a = 1,2, z = 1,2,3. (2.2) 



2 



Using the conventions of Ref. [|2T] we generically denote the couphngs of the neutral 
bosons to the fermions in the following by 



9^f = -Qf^ 9zf = ~—Qf^ 9zf = -—Qf + 



_ Mz _ -i 
gzzB = , 5'ZxH = , (2.3) 



where Q/ is the relative charge of the fermion /, and f = ±1/2 the third component 
of the weak isospin of the left-handed part of the fermion field /. The sine and cosine of 
the weak mixing angle are fixed by 



4-i-4-§. (2.4) 

In order to express the amplitude for the bremsstrahlung process e"'"e~ ttH7 in 
a compact way (see Section l^.3.ip . we use the Weyl-van der Waerden spinor technique 
as formulated in Ref. [^]. The helicity states for the massless electrons and massive 
top quarks are constructed as follows. Those for the massless incoming electron with 
momentum 

Pi = Pi (1, cos sin sin 0p J sin cos 6'p J (2.5) 
can be constructed from the Weyl spinor 

PM=V^r . Z " ' (2-6) 

V sm^ / 
and the corresponding momentum matrix reads 

Pi,AB = Pi,aP^,b- (2.7) 

The spinor p2,A and the momentum matrix P2 j^^ of the incoming positron are defined 
analogously. 

The spinors pi,A corresponding to the top-quark momentum 

fcf = {k^, |ki| cos0fc, sinOk^, |ki| sin0fc^ sin6'fci, |ki| cos6'fcJ (2.8) 

are defined as 

= i"'^'^: z ^ ] ' p^'^ = ^^^\ i ^ V (2-9) 

^ V sm / ^ V -e+"^'=i cos -f- J 

They obey the relation 

(P2P1) = mt, (2.10) 
where the spinor product is defined by 

= ^^e^'^^B = 01^2 - 02V'i- (2.11) 
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Figure 1: Tree- level diagrams for e"'"e ttH 

The spinors p'^ ^ corresponding to /c2 are constructed in the same way. The momentum 
matrices Ki and K2 corresponding the outgoing momenta ki and /c2 are decomposed into 
the spinors pi^A and p'^^^, respectively, 

^i,Ab = E Pi,APi,B, K^,AB = E Pi,AP'i,B- (2.12) 

i=l,2 i=l,2 

The Dirac spinors of the outgoing (anti-)top quarks are generically given by 

Ut{h) = {v\U, ^t(fc2)=(^^^) (2.13) 
with the actual insertions 

K.-rt-ji^'-T' T'^^' I '"'"f'^' T'^: (2.14) 

I (P2, Pi) for ri = -, [ (-P2, Pi) for T2 = +, 

i.e. we have (^77) = {rj'C,') = rrit. 

Although the Higgs boson has spin 0, it is convenient to express its momentum in 
terms of Weyl spinors (pi^A- These are obtained from the momentum analogously to 
(|2.9p . Then the momentum matrix of the Higgs boson can be decomposed as 

= E 'P^,A^^,B■ (2-15) 
i=l,2 

It is instructive to give the lowest-order amplitude for e"'"e~ —>■ ttH in this formalism. 
It receives contributions of the four diagrams shown in Figure ^ 

A^--i-2 ^ ^ _MX,^r,r, ^ j2 (^Mo''''^'^' + Mo''''^'^') , (2.16) 

y=7,z 

where the upper indices Z,x, t,t indicate the particle from which the Higgs boson is 
emitted. In the 't Hooft-Feynman gauge, the contributions corresponding to these four 
diagrams read 



_^z,.nr, ^ 2 ^7zc Pz{pi + P2) Pz{ki + k^) A^^^^^ (pi, ^i, k^), 
-Mr"" = -2e3^7L^?ZxH^xt^z(pi+P2)Pz(A:i + fc2)AJ.,.,(pi,P2,fci,fc2), 

^ _2e3^^^^HtPv(pi+P2)Pt(fci + A;3)Ar:,.,(pi,P2,fci,A;2), 
_yt<vW2 ^ -1e'g^y^g^,P^{vi^V2)P,{k2^k^)Al\^^{Vx.V2.k^M) (2-17) 
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with the auxihary functions 

^+,ri,r2(Pl>P2,/^l,fc2) = 9tt{p20* {PlO + 9zt{P2VT (PlV) , 

Al'^^^^^^{p^,P2,h,h) = gvt{P2vT{'^MPiV) + (e^3Pi)) 

+ 9^t{PiO{^MP20* - {P2K3V)[ 

AZ„rM^P2,h,h) = gUP20*{^rn,{p,e) + iv'KsPi)] 

+ 9vt{PiV){'^mt{p2V')* - {P2K3C) 

^-a,ri,r2(Pl>P2, A;i,fc2) = A^^^^^ (^2, Pi , , ^2) • (2.18) 

Here we used the shorthand notation 

{^|;'K^^|J) = ^'j^xt^'i^B = E (^>i)*(^0.) (2-19) 

i=l,2 

and abbreviated the propagators as 

Ptip) = ^^, PM = -r^. P,ip) = \- (2.20) 

p"^ — p"^ — M2 p 

Note that there is no need to introduce a finite width since none of the internal fines 
can become resonant for tfie pfiysical top-quark and Higgs-boson masses. Finafiy, tfie 
lowest-order cross section reads 

^o = YsJ ^"^3 ^(1 + 2i^-^)(l - 2P+a)|A<^|2, (2.21) 



where 

\M-o\' = N^ E E l-^r^l' (2.22) 

ri=±ir2=±i 

with the colour factor A^^'^ = 3 of the top quark, P± are the degrees of polarization of the 
e^ beams, and the integral over the three-particle phase space is defined by 

/ = (A / jSkf) (2-)'4p: + p. - g . (2.23) 

2.2 Virtual corrections 

2.2.1 Survey of one-loop diagrams 

The virtual corrections receive contributions from self-energy, vertex, box, and pen- 
tagon diagrams. The structural diagrams containing the generic contributions of vertex 
functions are summarized in Figure |21 The full set of pentagon diagrams is shown in 
Figure 01 
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eWteWHeWHe]^te]^t 



Figure 3: Pentagon diagrams for e"'"e — ^ ttH 



e l^Z t e t e W t 




Figure 4: e"'"e tt box diagrams 

The box diagrams for the e"'"e~tt vertex function are shown in Figure |3J the e"'"e^HH 
box diagrams in Figure and the e"'"e~Hx box diagrams in Figure IHl Figure [7| hsts the 
diagrams for the 7ttH and ZttH vertex functions. In Figure [71 as well as in Figure ITUI 
diagrams that can be obtained by reversing the charge flow of all charged particles are not 
shown. The e"'"e~ZH box diagrams can be found in Ref. [I2S1; and the e"'"e~7H diagrams 
can be obtained from a subset of those. 

The diagrams contributing to the 7tt and Ztt vertex functions are depicted in Figure |Hl 
the Htt and xtt diagrams in Figure IHl and Figure ^1 shows the jB.x and ZHx diagrams. 
The diagrams for the 7ZH, ZZH, and e"^e~H vertex functions are shown in Ref. [1^. 
Those for the 77H vertex function can be obtained from a subset of the 7ZH vertex 
diagrams. The e"'"e~7 and e"'"e~Z diagrams and most of the self-energy diagrams can be 
found in Ref. [ |21] . The diagrams for the top-quark self-energy are listed in Figure ^2 The 
ZH-mixing energy and the 7HH and ZHH vertex functions vanish owing to CP symmetry. 

All pentagon and box diagrams are ultraviolet (UV) finite, as well as the 77H vertex 
function. Since we neglect the electron mass whenever possible, also the e"'"e~H vertex 
function is UV finite. The remaining vertex functions and self-energies are UV divergent, 
and the corresponding counterterm diagrams have to be included. 
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Figure 5: e"'"e HH box diagrams 



e W H e W X 




e \Y X e \Y H 



Figure 6: e"'"e Hx box diagrams 

2.2.2 Calculational framework 

The actual calculation of the one-loop diagrams closely follows the strategy of Ref . [ 122] , 
where the 0{a) corrections to e+e~ uuB. have been calculated, i.e. it has been car- 
ried out in the 't Hooft-Feynman gauge using standard techniques. The Feynman graphs 
have been generated with FeynArts [|2H1 and are evaluated in two completely independent 
ways, leading to two independent computer codes. The results of the two codes are in 
good numerical agreement (i.e. within about 12 digits for non-exceptional phase-space 
points). Apart from the 5-point integrals, the tensor coefficients of the one-loop inte- 
grals are algebraically reduced to scalar integrals with the Passarino-Veltman algorithm 
[ 1^ at the numerical level. For the evaluation of the tensor 5-point functions the direct 
reduction to 4-point integrals of Ref. [|2Z1 has been used. This method avoids leading 
inverse Gram determinants which cause numerical instabilities at the phase-space bound- 
aries. The scalar integrals are evaluated using the methods and results of Refs. [ |2Il I2H] , 
where UV divergences are regulated dimensionally and infrared (IR) divergences with an 
infinitesimal photon mass m^. The soft and coUinear singularities of some scalar 5-point 
integrals have been checked against the general result of Ref. [ 121] . The renormalization 
is carried out in the on-shell renormalization scheme, as e.g. described in Ref. [I21j. 

In the first calculation, the Feynman graphs are generated with FeynArts version 1.0 
[ I25j. With the help of Mathematica routines the amplitudes are expressed in terms 
of standard matrix elements (SME) and coefficients of tensor integrals, as described in 
the appendices in more detail. The output is processed into a Fortran program for the 
numerical evaluation. The second calculation has been made using FeynArts version 3 
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Figure 8: Diagrams for 7tt and Ztt vertex functions 




Figure 9: Diagrams for Htt and xtt vertex functions 

[ Eni for the generation and FormCalc [ |^ for the evaluation of the amplitudes. The 
analytical results of FormCalc in terms of SME and their coefficients were translated into 
C++ code for numerical evaluation. 

As a check of gauge invariance the calculation of the virtual corrections has been 
repeated using the background-field method [ |321, where the individual contributions 
from self-energy, vertex, box, and pentagon corrections differ from their counterparts in 
the conventional formalism. The total one-loop corrections of the conventional and of the 
background-field approach were found to be in perfect numerical agreement. 

Finally, the contribution of the virtual corrections to the cross section is given by 

Sa^,rt = ^^ /d$3 E \{l + 2P-a){l-2P^a)N^ ^ ^ 2 Re {A^^- (A^r^D , 

CT=±2 Ti=±^T2=±2 

(2.24) 

where JK4Y^^ denotes the one-loop helicity amplitudes. 
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Figure 10: Diagrams for 7HX and ZH^ vertex functions (Diagrams with reversed charge 
flow in the loop are not shown.) 




t b t b 



Figure 11: Diagrams contributing to the top-quark self-energy 

2.3 Real photonic corrections 

2.3.1 Matrix- element calculation 

The real photonic corrections are induced by the process 

e-(pi, ai) + e+(p2, ^2) — > t{ki, n) + t{k2, T2) + H(/c3) + -f{k, A), (2.25) 

where k and A denote the photon momentum and helicity, respectively. The Feynman 
diagrams of this process are shown in Figure IT^ We have evaluated the corresponding 
helicity matrix elements Ai'^'^'^'^^''^ using the Weyl-van der Waerden spinor technique as 
formulated in Ref . [ . The amplitudes for the helicity channels with cji = (J2 vanish for 
massless electrons, and we define a = ai = —a2 as above. The treatment of IR (soft and 
collinear) singularities, where the finite electron mass is needed, is described below. 
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Figure 12: Feynman diagrams for e^e ttH7 



The amplitude for e+e —>■ ttH^y can be separated into 12 contributions that are 
invariant with respect to gauge transformations of the external photon, 

j^aTiT2X _ j^Ze,aTiT2X _|_ j^Zt,aTiT2X _|_ ^^xe,o'nT2A _|_ ^^xt:0"m"2A 

_j_ ^J^Vte,o-TiT2\ _j_ j\^Vtt,crTiT2X _j_ j\^Vtc,aTiT2X _|_ j\^Vtt,(TTiT2X^ ^2 26) 



y=7,z 



The last upper particle index in etc., indicates to which charged fermion line (e or t) 
the photon is attached. The diagrams in the first line of Figure IT^ belong to Ai'^^ and 
TW^*, those in the second line to A4^^ and A4^^, those in the third line to A4Y^^ and A4 



and those in the last line to and 

In the 't Hooft-Feynman gauge, the individual contributions read 

-V2e^ Qt ^?ze ^?zzH Pz(pi + P2) Pzih + k2 + k) 
X All^^^^{pi,p2,ki,k2,k), 



Vtt 
7 ' 



Zt,crrir2 A 
7 



(2.27) 
(2.28) 
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j^xe,ar,r,x ^ 2^2 c^g^^^ gz^B 9xt PziPl +P2-k) Pz{h + k2) 

xAf^^^^^{p,,P2,h,k2,k), (2.29) 

^Xt,'rr,r,X ^_2^e^ Q^g-^ ^^xH Q^t Pz (Pl + P2) Pz {h + ^ + k) 

X Af^^^^{p^,P2MM.k), (2.30) 

^Vte,ar,r,X ^ ^^e^g^^^g^, Py{p, + P2 - k) Pt{ki + ^3) 

xAl^^^^^{p,,P2,k,,k2,k), (2.31) 

^vtt,.r,r,x ^ -2V2e'Q,g^y^gBtP^,{pi+P2)Alll^^{p,,p2MM.k), (2.32) 

^V^.,ar,r,X ^ p^^^^ + P2 - k) ^(^2 + k^) 

xAl~^'^^^^{p,,P2,k,,k2,k), (2.33) 

^v-u,ar,r.x ^ -2V2 Q, g^^ gut Py(pi + P2) A^^l^^(p,, p2, k,, k2, k) (2.34) 
with the auxihary functions 

A A. ^(pi,P2, ki,k2,k)^ -9zt{Pi0^j — VTi — ^ 9zt{PiV) 



{kpi){kp2) ^^tu-^'/z {kpi){kp2) 

{kK2Kik) 



Alr,,r,,M^P2,kl,k2,k) = {gl{p20* {PlO + Qz^T M) ^^^^ _ 



+ {kp2)' 



( gtMY 

\ {ki ■ k) 



{piO - 



{k2 ■ k) 



Al%^^,^^^{p„P2,k„k2,k) = - ii'v)) 



{PiV)j , 

{p,{P2-K)K3Pir 
{kpi){kp2) 



Af^r^^r^^+{pi,P2,ki,k2,k) = {P2K3P1) 



AVtC 



4:{ki-k){k2 
iPi,P2, ki, k2, k) = 



1 1 

+ 



1 



{kpi){kp2) 

+ g+,{p,e){2m,{C{P2 - K)p,) + (^^(P^ - K)K3r^Y) 



2{ki ■ k) 2{k2 ■ k) 
gvt{v'{P2 - K)p,){2m,{pirj) + {iK^pi)) 



Ay;;^^^^^^{pi,P2,ki,k2,k) 



{kKsK,k)Pt{ki + ks + k) - 



{kK2Kik) 



X 



2{k2 ■ k) 

gUPiO (2mt(p2e)* - {P2Ksv)) + 9vt{P2vT C^rritipiv) + i^KsPi)) 

{kO 



2{ki ■ k) 



Pt{ki + k3 + k) 2g+,mt(j>iO{kp2Y - 9vx{P2vT {kK^Pi) 
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A 



Vie 



- {kp2r{piOPt{ki + h) Ptih + h + k) g+, [2m,{kO* - (kK^v)] 

1 



Xpi,P2, ki, k2, k) 



{kpi){kp2) 

+ gy,{p,r]){2m,{r]'{P2 - K)p,) + {p.iP, - K)K,ir 



gUm - K)p,)(2m,{prO + {r^'K.pr) 



2{k2 ■ k) 



Pt{k2 + ks 



{kKsK2k)Pt{k2 + ks + k) + 



{kK2Kik) 



+ 



2{ki ■ k) 

9vt{P20* (2mt(pi^') + WK3P1)) + gvtiPiv) {2mt{p2v'y - {P2K3O) 
{kvT 



2{k2 ■ k) 
{kO*{kp2) 



Pt{k2 + ks + k) [2gy,mt{pir]){kp2)* - g^t{P2iY {kK^Pi) 
Pt{k2 + A;3) g^, [2m,{p,0 + (V^aPi)] 



2{ki ■ k) 

+ {kp2)*{piv)Ptik2 + ks) Pt{k2 + k3 + k) gy, [2m,{kr]'y - {kKsO] (2.35) 
and the symmetry relations 

A-a,TuT2,xiPi^P2,ki,k2,k) = -A-^^^^^{p2,pi,ki,k2,k), 
A^L^^^^^^_^{pi,P2,ki,k2,k) = sgn{TiT2)A^J^^^^{p2,pi,k2,ki,k)*, B = Z,x, / = e,t, 
Al^Ir2,-Ti,-xiPuP2, ki, k2, k) = sgn{TiT2)Al^l^^^{p2,pi, k2, ki, k)*, / = e,t, 



iPi,P2, ki, k2, k) = sgn{TiT2)Al^(^ Jp2,pi, k2, ki, k)\ / = e,t 



(2.36) 

Note that interchanging ki ^ k2 means also interchanging r) ■ir^ ^' , ^ ^ r]' . Besides (I2.19j) 
we used the following auxiliary quantities in the previous definitions: 

{kK2K,k) = k^Ki''K,,,^k^= y: {kp[r{p,p[){p,kr, 

{kK,K,k) = kj,Ki^K,j,ck'^ = {kcl>i)*{p^ct>,){p,k)\ 
{kK^K2k) = kj,Ki^K,j,ck^ = Y {kcl>i)*{p'^ct>,){p'^k)\ 

= E ((V^'P2)*(0,P2)(0,^)* - (V^'fc)*(0,fc)(0,V^)*) . (2.37) 
i=i,2 

The Weyl spinor k^ and the momentum matrix Kj^^ of the photon are defined from the 
photon momentum k analogously to those of the massless electron in ()2.fjj) and ()2.7j) . 
The contribution of the radiative process to the cross section is given by 



a\2 

7I 



(2.38) 
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with 



\M^,\' = A^t E E E l-^r"'!'' (2-39) 

A=±l ^,=±1 ^2=±i 

and the phase-space integral 

I^^^-Ii^^ {Ujwm) (2.40) 

The spin- averaged squared matrix element obtained from Eqs. ()2.26p - ()2.36p has been 
successfully checked against the result obtained with the package Madgraph [EH] numer- 
ically. 

2.3.2 Treatment of soft and collinear singularities 

Without soft and collinear regulators the phase-space integral (j2.38j) diverges in the 
soft (fco 0) and collinear [pak —>■ 0) phase-space regions. In the following we describe two 
procedures of treating soft and collinear photon emission: one is based on a subtraction 
method, the other on phase-space shcing. In both cases soft and collinear singularities 
are regularized by an infinitesimal photon mass and a small electron mass, respectively. 

Comparing both methods we find numerical agreement for total cross sections and 
distributions within integration errors. The results presented in this paper were obtained 
with the subtraction method, and the statistical error is typically 0.1% or below for the 
total cross section. 

(i) The dipole subtraction approach 

The idea of so-called subtraction methods is to subtract a simple auxiliary function 
from the singular integrand of the bremsstrahlung integral and to add this contribution 
back again after partial analytic integration. This auxiliary function, denoted |A^subP 
in the following, has to be chosen in such a way that it cancels all singularities of the 
original integrand, which is jA^^P in our case, so that the phase-space integration of the 
difference can be performed numerically, even over the singular regions of the original 
integrand. In this difference, Ai-y can be evaluated without regulators for soft or collinear 
singularities, i.e. we can make use of the results of the previous section. The auxiliary 
function has to be simple enough so that it can be integrated over the singular regions 
analytically, when the subtracted contribution is added again. This part contains the 
singular contributions and requires regulators, i.e. photon and electron masses have to be 
reintroduced there. Specifically, we have applied the dipole subtraction formalism, which 
is a process-independent approach that was first proposed [,34^ within QCD for massless 
unpolarized partons and subsequently generahzed to photon radiation of massive polarized 
fermions in Ref. [EH]- In order to keep the description of the method transparent, we 
describe only the basic structure of the individual terms explicitly and refer to Ref. [EH] 
for details. 

In the dipole subtraction formalism the subtraction functions are constructed from 
contributions related to ordered pairs of charged fermions. These two fermions are called 
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emitter and spectator, respectively, since by construction only the kinematics of the emitter 
leads to coUinear singularities. In the following, we consider only the case for unpolar- 
ized final-state fermions. The subtraction function receives four kinds of contributions 
depending whether the emitter (first index) and spectator (second index) are chosen from 
initial (labelled by a,b) or final state (denoted by 

l-Msubr= E \MU,at\'+ E {\MUJ' + \MUJ')+ E l-^sV./ (2-41) 



a, 6=1, 2 



a=l,2 
1=1,2 



with the contributions 

,ia I 

l-'^sub,jjl ~ 



2 (sub) 



:-l)"+*+^Qte'^7iff (Pa, k) MUpi,P2, ku k2, k,) 

i-ir+'+'Q,e'gt''\h,Pa. k) \MUpuP2. h, h] 

Qle^ 9fj'^\ki, kj, k) Mo{pi,p2, h, h, h] 
The unpolarized dipole functions for final-state emitters read 

(sub) , 



(2.42) 



9ia,+ {ki,Pa,k) + Qia,- {ki,Pa,k), 



9ia 

9ij kj, k) = g\j^^'{ki, kj, k) + g\j _'{ki, kj, k). 



(sub) , 



(2.43) 
Note 



The dipole functions g^^^^^ are defined in Eqs. (3.22), (A.l), and (4.4) of Ref. 
that the spin-flip dipole functions g^^^ and ga^^^ vanish in the case of small fermion 
masses, i.e. rrie — * 0. The arguments pa, ki of the tree-level matrix elements Mq in 
subtraction functions depend on the momenta Pa, ki, k. The mapping from Pa, ki, k into 
Pa, ki is defined by Eqs. (3.25) and (3.26) of Ref. [ESI for -^sub,ab' by Eq. (3.12) for M^^^^^i 
and M^^^^ia, and by Eq. (4.5) for M^^^^ij. 

The subtracted contribution can be integrated over the (singular) photonic degrees 
of freedom up to a remaining convolution over x (= Xab,Xia)- In this integration the 
regulators and mg must be retained, and the soft and collinear singularities appear 
as logarithms in these mass regulators. This non-trivial step has, however, to be done 
only once and for all, and the needed results can be found in Ref. [ESI- The integrated 
subtraction function reads 



asnhiPl,P2,P-,P^ 



a 
2^ 



d^E E 



.T=± a=l,2 



E^, 

6=1,2 
bj^a 



(sub) 
ab,T ' 



+ E(-i 



ya+i+l 



X] 



i=l,2 



+ E 



a=l,2 
i = l,2 



+ 



a 
2^ 



E E 

.T=± a=1.2 



dcT, 



(a,r). 



l)"+^ + 'gt [Q; 

E^(sub) 



(sub) 



{tai 1 x) 

+ E(-i)"^^^'QtGi?^(t. 



6=1,2 
6/a 



j=l,2 



+ / dao 



(sub) , 



a = l,2 
1=1,2 



i,J = l,2 



(2.44) 
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with the modified lowest-order cross sections 

da^^'''\x) = dao{xpi,p2,TP_, P+), da^^'''\x) = dao{pi,xp2, P-,tP+) (2.45) 
and the unpolarized distributions Q'^^^^^ and end-point contributions G^'^^^^ 

Gia {taiiX) = Gia^+ {taiiX) + ^jq, _ (tai,x), 
^(sub)/, \ _ ^(sub)/, \ _|_ (<^(sub) / , X 

Gtf^{s.,) = G^{s.,) + Gt-^{s.,). (2.46) 

In our case, the relevant definitions for Q^^^^^ and G^^^^^ are given in Eqs. (3.32), (3.33), 
(A. 2), (A. 4), and (4.10) of Ref. The [. . .]+ prescription is defined as usual. 



dx 







fix) q{x)= / dxf{x)[g{x)-g{l)]. (2.47) 







In summary, the phase-space integral ()2.38|1 in the dipole subtraction approach reads 

^ / E ^(1 + 2^-^)(l - 2i^+^) - l-^^ubP] + ^sub. (2.48) 



(ii) The phase-space-slicing approach 

The idea of the phase-space-slicing method is to divide the bremsstrahlung phase space 
into singular and non-singular regions, then to evaluate the singular regions analytically 
and to perform an explicit cancellation of the arising soft and coUinear singularities against 
their counterparts in the virtual corrections. The finite remainder can be evaluated by 
Monte Carlo techniques. For the actual implementation of this well-known procedure we 
closely follow the approach of Ref. [ EE] • We divide the four-particle phase space into soft 
and collinear regions by introducing the cut-off parameters AE and A6, respectively, and 
decompose the real corrections as 

day = ddsoft + ddcoii + day 

,finite- 

(2.49) 

Here dcxsoft describes the contribution of the soft photons, i.e. of photons with energies 
fco < AE in the CM frame, and dacoii describes real photon radiation outside the soft- 
photon region (fco > AE) but collinear to the e^ beams. The collinear region consists of 
the two disjoint parts < 9y < AO and tt — AO < Oy < n, where Oy is the polar angle of 
the emitted photon in the CM frame. The remaining part, which is free of singularities, 
is denoted by da^^gnite- 

In the soft and collinear regions, the squared matrix element |A^^^|^ factorizes into the 
leading-order squared matrix element |A^oP ^ collinear factor. Also the four- 

particle phase space factorizes into a three-particle phase space and a soft or collinear 
part, so that the integration over the singular part of the photon phase space can be 
performed analytically. 

In the soft-photon region, we apply the soft-photon approximation to |A^J^p, i.e. the 
photon four-momentum k is omitted everywhere but in the IR-singular propagators. In 
this region da^ can be written as [ EI] 

a f d^k ( pi p2 ^ ki ^ ^2 ^ ^ 
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This can be decomposed into an ISR part, a final-state-radiation (FSR) part (proportional 
to Qt) and the ISR-FSR-interference part (proportional to Qt), 



do". 



soft 



do". 



soft, ISR 



+ da, 



soft,FSR 



softJnt- 



(2.51) 



The explicit expressions for the soft-photon integrals can be found in Refs. [ |^ |2H1 • 
For our purpose it is sufficient to keep the electron mass only as regulator for the coUinear 
singularities. In this limit we obtain for the ISR part 



da. 



softJSR 



-dan 



a 



TT 



2 In 



/2AE\ 



In 



In 



1 



In^ 



nit 



The FSR part reads 



a 



ddsoft.FSR = -dao—i^Qt [ht + Iti - 2/i 



47r 



tt I 



(2.52) 



(2.53) 



with the massive integrals la and lij from Refs. [ El I2H] • Finally, the interference contri- 
bution is given by 



da. 



soft,int 



-dan 



a 
'2^ 



■Qt [/c 



-^c-t - -^e+t + It 



+t 



(2.54) 



in terms of the soft integral with a light particle a and a massive particle j, 



■■aj 



27r 



In 



f 2AE \ ^^f{2paPjf 

Li2 (1 



n 2 

Pa 



--In^ 



{P°j + \Pj\)PaPj, 



Li2 



(P,° + |P.I)V 4 

Pa{p' + \pj\y 



~6 



PaPj 



(2.55) 



In the collinear region, we consider an incoming e^ with momentum Pa being split into 
a collinear photon and an e"^ with the resulting momentum xpa after photon radiation. 
The cross section for hard photon radiation (/eq > AE) in this region reads 



acoi\ipi,P2,P-,P^ 



a 
2tx Jo 



l-2A£;/v^ 



do; E 

r=± 



(coll), 



S,X] 



da^^^^\x) + 



da, 



X] 



where fl2.45|) is used, and 



X 



In 



4m? 



,(coll) 



s,x) 



X. 



(2.56) 



(2.57) 



Finally, the finite cross section da^^finite is defined by imposing cuts on the bremsstrah- 
lung phase space ()2.4m) . i.e. a photon-energy cut, > AE, and a cut on the angles 
between the photon and the beams, < < vr — A6. 

For the numerical evaluation the values AE/y/s = 10^^ and A^^ = 10^^ have been 
chosen. For ^/s = 500 GeV and Mh = 140 GeV as well as for = 1 TeV and Mh = 
115 GeV, we have checked that the cross section changes only within integration errors 
when varying AE within 10~^ GeV < AE < 10~^ GeV for fixed A6 = 10~^ and when 
varying AO within 10~^ < AO < IQ-^ for fixed AE = IQ-^ GeV. 
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2.4 QCD corrections 

The 0{as) QCD corrections can be obtained from the photonic FSR part of the 0{a) 
correction, i.e. the part proportional to Q^, by the replacement 

Qla ^ Cpa^ifi^) = ^as(/i'), (2.58) 

where Qt = 2/3 is the relative charge of the top quark. The real part of the FSR 
corrections results from diagrams where a photon is emitted from top quarks. The cor- 
responding matrix elements Ai'^-^ are given in Eqs. ()2.28|) . ()2.30p . ()2.32p . and ()2.34|) . 
Following Ref. [112|, the QCD renormalization scale n is set to the CM energy ^/s, and 
the running of the strong coupling is evaluated at the two-loop level (MS scheme) with 
five active flavours, normalized by as{M^) as given in 1)3.11) . For = 500 GeV, 800 GeV, 
and 1000 GeV the resulting values for the strong coupling are given by as(M|) = 0.09349, 
0.08857, and 0.08642, respectively. 

2.5 Initial-state radiation beyond 0{a.) 

The emission of photons coUinear to the incoming electrons or positrons leads to cor- 
rections that are enhanced by large logarithms of the form ln{ml/s). In order to achieve 
an accuracy at the few 0.1% level, the corresponding higher-order contributions, i.e. con- 
tributions beyond 0{a), must be taken into account. These are included in our calculation 
through a convolution of the leading-order cross section with leading-logarithmic structure 
functions precisely in the same way as described in Ref. 

2.6 Monte Carlo integration 

The phase-space integration is performed with Monte Carlo techniques in both com- 
puter codes. The first code employs a multi-channel Monte Carlo generator similar to 
the one implemented in RacoonWW [EH 1^^] and Lusifer [EO], the second one uses the 
adaptive multi-dimensional integration program VEGAS [HT]. 

3 Numerical results 



3.1 Input parameters 

For the numerical evaluation we use the following set of SM parameters [ |42j : 



Gn 


= 1.16639 X 10-^GeV-^ 


a(0) 


= 1/137.03599976, 


«s(Mz) 


= 0.1172, 


Mw 


= 80.423 GeV, 


Mz 


= 91.1876 GeV, 






rUe 


= 0.510998902 MeV, 




= 105.658357 MeV, 


rrir 


= 1.77699 GeV 




= 66 MeV, 


rric 


= 1.2 GeV, 


mt 


= 174.3 GeV, 


rrid 


= 66 MeV, 


rus 


= 150 MeV, 




= 4.3 GeV, 



(3.1) 

which coincides with the one used in Refs. [1201 123- Since we employ the on-shell renor- 
malization scheme, the weak mixing angle is fixed by ()2.4)1 . Because quark mixing is 
irrelevant for the considered process we use a unit quark-mixing matrix. 
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We evaluate amplitudes in the so-called G^-scheme, i.e. we derive the electromagnetic 
coupling a = e^/(47r) from the Fermi constant according to 

TT 

This procedure, in particular, absorbs all sizable mass effects of light fermions other than 
electrons in the coupling a^^, and the results are practically independent of the masses of 
the light quarks. The masses of the light quarks are adjusted to reproduce the hadronic 
contribution to the photonic vacuum polarization of Ref. [ESI- In the relative radiative 
corrections, we use a(0) as coupling parameter, which is the correct effective coupling 
for real photon emission. We do not calculate the W-boson mass from but use its 
experimental value as input. 

3.2 Results on total cross sections 

Numerical results for total cross sections including individual contributions to the ra- 
diative corrections have already been presented in Ref. [ 120] • The contributions discussed 
there were the photonic corrections originating from virtual photon exchange and real 
photon radiation, the beyond 0{a) or higher-order ISR corrections that are included us- 
ing the structure-function approach as described above, and the remaining non-photonic 
weak corrections. All these have been combined into the complete electroweak corrections. 
Furthermore, the QCD corrections have been considered. 

Here, we complement this by a splitting of the weak corrections into the purely 
fermionic loops together with the fermionic part of the counterterms and the remain- 
ing weak bosonic contributions. Instead of the separate photonic and higher- order ISR 
contributions considered in Ref. [120]; we only look at their sum, denoted by QED. The 
corresponding relative corrections are shown as a function of the CM energy in Figure ^1 
for Higgs-boson masses of Mh = 115 GeV and 150 GeV. 

The fermionic corrections are about +10% and depend only weakly on the CM energy. 
While the weak bosonic contribution is also around +10% close to threshold, it falls off 
rapidly with increasing CM energies, eventually reaching about —20% at an energy of 
1.5 TeV. Therefore, for energies above 600 GeV, the fermionic and weak bosonic con- 
tributions partially cancel when combined to form the complete weak corrections. The 
QED corrections are about —40% at threshold rising to a few per cent at 1.5 TeV. In 
the electroweak corrections, both QED and weak contributions partially compensate each 
other, yielding relative corrections of about —20% at threshold, and —3% and —7% for 
Mh = 115 GeV and 150 GeV, respectively, at 1.5 TeV. 

3.3 Results on distributions 

The distributions in the Higgs-boson energy i?H and polar angle are shown in 
Figure HH for Mh = 115 GeV and ^/s = 800 GeV, where the cross section is nearly 
maximal for small Higgs-boson masses. All energies and angles used in the distributions 
here and below are defined in the CM frame. The curves denoted with electroweak contain 
only the electroweak corrections and (in case of the absolute distributions) the lowest-order 
contribution, whereas in the curve labelled full the QCD corrections are included as well. 
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Figure 13: Relative corrections for Mh = 115 GeV (l.h.s) and Mh = 150 GeV (r.h.s) 



The distribution in the Higgs-boson energy is a broad continuum whose maximum around 
190 GeV gets slightly shifted down to smaller energies by the electroweak corrections but 
is left unchanged by the QCD corrections. At the upper end of the £'h spectrum the tt 
system becomes non-relativistic, and the QCD correction is dominated by the Coulomb 
singularity [ [T21 CH HE] induced by soft-gluon exchange in the tt system; the rising QCD 
correction with increasing i?H reflects the transition to this singular region. Although 
photonic FSR shows the same behaviour, this effect is totally overruled by ISR in the 
QED corrections. For increasing £^h these corrections become more and more negative 
because of the reduced phase space after ISR, which prohibits a cancellation of large 
IR-sensitive virtual corrections by the real corrections. For the chosen CM energy and 
Higgs-boson mass, the purely weak corrections depend only weakly on the Higgs-boson 
energy i?H and are only at the per-cent level. Note, however, that the smallness of the 
corrections for the chosen setup is mainly due to the transition from larger positive to 
negative corrections with increasing CM energy, as can be seen in Figure IT^ Thus, for 
lower and higher CM energies sizeable weak corrections to the £'h spectrum arise. 

The distribution in cos ^h reaches maxima in the forward and backward directions at 
leading order. The inclusion of the radiative corrections does not sizeably distort this 
shape, and essentially only the normalization is changed. However, viewed separately the 
weak and photonic contributions show a small dependence on ^h which cancels in their 
sum. 

Figure ^1 shows the distributions in the top-quark energy i?t and in the cosine of the 
angle between the top quark and the incoming electron, cos6't, for = 800 GeV and 
Mh = 115 GeV. The Ex, distribution has a broad peak at high energies, which is smeared 
out by the QED and QCD corrections towards lower energies. While the weak corrections 
vary only slowly with E't, QCD and QED corrections decrease rapidly and reach large 
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Figure 14: Distribution in the Higgs-boson energy En (upper panels) and the Higgs-boson 
polar angle 6'h (lower panels) and corresponding relative corrections for i/i = 800 GeV 
and Mh = 115 GeV 
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Figure 15: Distribution in the top-quark energy (upper panels) and the top-quark 
polar angle 9t (lower panels) and corresponding relative corrections for ^/s = 800 GeV 
and Mh = 115 GeV 
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Figure 16: Distribution in the photon/gluon energy E^/g and in the photon/gluon polar 
angle cos in the radiative processes e+e" — ttH + '~f/g for ^/s = 800 GeV and Mh = 
115 GeV 



negative values for large Et. This is due to the fact that for large Et, the available 
phase space for the emission of a real photon or gluon is restricted and the large negative 
virtual corrections can only be partially compensated. The distribution in the top-quark 
production angle is maximal in the forward (e~) direction. The relative QCD corrections 
show only a small angular variation, but the weak and QED contributions, and hence the 
electroweak corrections, depend on cos6't, weakening the maximum somewhat. Owing to 
CP-symmetry the energy distribution of the anti-top quark is identical to the top-quark 
distribution. The angular distribution for the anti-top quark can be obtained from the 
one of the top quark by the replacement cos6't — > — cos6't. 

Finally, the spectra of the photon energy E^ and of the photon polar angle cos 6^ of the 
radiative process e+e~ —>■ ttH-l-7 are shown in Figure ITHl together with the corresponding 
distributions for the gluon in the process e"'"e~ ttH -|- g, again for = 800 GeV and 
Mh = 115 GeV. In order to render the photon/gluon visible, we impose angular and 
energy cuts of 

9{-f /g, beam) > 5°, E^/g > 5 GeV. (3.3) 

The photon and gluon energy spectra are dominated by the IR pole in the region of small 
energies. Owing to the coUinear singularity of ISR, the photon angular spectrum shows 
a strong peaking behaviour in the forward and backward directions. In contrast, for the 
chosen set-up the gluon radiation does not show a dependence on the gluon production 
angle. The reduction in the first and last bins is due to the angular cut. 
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4 Summary 



If the Higgs boson is not too heavy, an experimental analysis of the process e"'"e^ ttH 
at a future high-energy e"'"e~ linear collider allows for a determination of the top-quark- 
Higgs-boson coupling at the level of 5% or better. At this level of accuracy, the inclusion 
of both QCD and electroweak radiative corrections is important. 

We have presented the calculation of the 0{a) and 0{as) radiative corrections to the 
process e"'"e~ — > ttH in the Standard Model in some detail. Particular attention has been 
paid to the treatment of the soft and collinear singularities in the real photonic (and 
gluonic) corrections. 

Numerical results on the corrections to total cross sections are discussed, complement- 
ing existing work in the literature; in particular, the interplay of weak fermionic, weak 
bosonic, and photonic corrections is investigated. We find that all these contributions 
are typically at the level of 10% and that the different contributions partially cancel each 
other. Moreover, we have presented numerical results for the electroweak corrections to 
differential cross sections, including energy and angular distributions of the final-state 
particles. In the G^-scheme the size of the corrections typically reaches the 10% level 
and even exceeds it near edges of phase space, underlining the importance of a thorough 
understanding of the radiative corrections. 

Appendix 

A Standard matrix elements 

For the calculation of the loop corrections it is convenient to separate the matrix ele- 
ments into scalar invariant coefficient functions containing the loop integrals and standard 
matrix elements (SME) containing all tensorial and spinorial objects and the dependence 
on the helicities of the external particles [ |^ . 

To introduce a compact notation for the SME, the tensors 

rJl^A..^.^,...} = Utih) {1, 7^, 7^7^, . . .} UJrVt{k2) (A.l) 

are defined with uj± = {1 ± 75)/2 and an obvious notation for the Dirac spinors Ue-{pi), 
etc. Furthermore, symbols like Tp are used as shorthand for the contraction F^p'^. In 
close analogy to the calculation [1^ of QCD corrections to the related process qq — > ttH, 
we introduce the following set of SME, 

,cr,fj, -pttjT 

{fi,fj.pi,fip2,fipiP2}^ 
,a,^k\k2 ptt.T 

{/^,/ipi,/ip2,AtpiP2}' 

,a,k\ pttjT 

{l,pi,P2,PlP2}' 

,(7,^2 pttjT 

{l,Pl,P2,PlP2}' 
,a,fiuki ptt,T 

{fJ,U,^lupi,^lup2,fJ.upiP2}^ 



•'^'{1,2,3,4} — ^ 

•'^'{5,6,7,8} — 

•^{9,10,11,12} = 

•^{13,14,15,16} = 

•'^'{17,18,19,20} ~ 
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»^(Tr -pe+e ,cr,fj,uk2 ptt,T 

{21,22,23,24} {fj,i>,^upi,^up2,fJ.i'pip2}'' 

i^ACTT _ T^e+e- ,a,tMup ptt,T ( A O') 

•'^l {25,26,27,28} — {iiup,tiuppi,tiupp2,tJ.uppip2}- l^-^J 

The four-dimensionality of space-time implies that the SME Ai'J'^ are not all indepen- 
dent; there are linear relations among them. A simple way to derive the relations with 
real coefficients is provided by the following trick. In four dimensions the metric tensor 
can be decomposed in terms of four independent orthonormal four- vectors rii, 

g^P = n'^,4-jZnrnl (A.3) 
1=1 

where Uk ■ rii = g^i. Two convenient choices (j = 1, 2) for the vectors rii are given by 



< = J r r ' 2 + %? + %2 

Y tijt2j - mis \ ^ s s ^ 

2 

a ,Cf/37'5„ n n Z af^'yS 

3 „ ^ - „ . 



n 



e"^^Vo,/3ni,^n2,5 = , _ _ ^ e°^^Vi,ffP2,-,fci,^ (A.4) 

s{tijt2j - mfs) 



with 

tij = tij - ml, i,j = 1, 2, (A.5) 

and e'^^')''' (^^0123 _ j^-^^ denoting the totally antisymmetric tensor. Inserting this decom- 
position for both j = 1,2 in all contractions between different Dirac chains according to 
-pe+e ,(T,o-ptt,T _ -pe+e gO'f^Y^^''^ , ctc, and Subsequently using the Dirac equation and the 
Chisholm identity 



le 



reduces all 112 SME M.'^'^ to 16. Further relations with complex coefficients result from 
the direct application of the Chisholm identity ()A.6|) to structures like ^l^kip" subse- 
quently using the decomposition 

4 

gP'^ = ^ fz-') , . 2pfpJ, Z,, = 2p,pj, PS = h, p, = h, (A.7) 

to separate all contractions between e tensors and Dirac chains via .'^^0/375 _ 

pe+e ,agpS^ai3-y^^ ^p2kipt ' particular, this leads to 



pe+e ,cT p2kipiS _ -pe+e ,u p5 p2k\Pi _ rfy 



p 

with 



•n^"^3+rr"(^"/44 



(A.8) 



X = ^v^P2k,k2 ^ e^^'"'pi,^p2,uKpk2,.. (A.9) 

In the inverse matrix {Z^^), the determinant det(Z) occurs, which can be identified with 
det(Z) = -16X2. 
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Altogether, the hnear relations reduce the set of SME to 8 SME. It is convenient to 
express all SME in terms of Ml^, Mt"^, and Mf^ (-^2^ and ^^3='' vanish), 

Mi+ = ttiMt^ + biMf- + CiMt^ + diMf, 

M-- = a*Mr + KM-, + + <Mr + rf*^3-+, 

Mt = eMt^ + fMt- + gMt^ + hMt~, 

M;+ = e*Mr + f:Ml^ + 9*Mr + h*Ms+. (A.IO) 

The non- vanishing coefficients a^, etc., for the dependent SME read 

_ D*ml - 2si2ii2i22 _ mt(B*tn - 2mlsti2) 

_ t22('2mlstii - Ctu) _ rrit _ hi 

~ 2E2 ' ~ T' ~ "T' 
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Ol9 — — ^ , O21 — — , 022 — — 4ri2, 

023 = -4^22, 025 = 16, (A-H) 
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2Eis ' " 2E2S ' 2s' 
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_ ill , _ B*tn-2mlsti2 , _ A*mt 
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£/l 

h27 = 16, (A.18) 
where the following shorthands are used, 

A = SS12 - tl2t21 - ^11^22 + 4iX, B = SS12 + ti2t21 - tilt22 + 4iX, 

C = SS12 - tuhl + iiit22 + 4iX, D = ssu + t^hi + ^11^22 + 4iX, 

Ei ^ -mls + iuhi, i = l,2, (A.19) 

and 

S12 = S12 - 2ml. (A.20) 
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Finally, we give the explicit expressions of the independent SME Mi^, Mf^, and 
Aif^ in terms of Weyl-van der Waerden spinor products as introduced in Section I^ITI 

Mt+ = 2(M)*(M'), Mr = 2(pir/')*(P2r/), 

Mt- = 2{p20*{P2Vr{PiP2), M^^ = -2{piP2r{p20{P2V), (A.21) 

where the insertions for the spinors ^^'^ and r/^'-* for the polarizations ri^2 are specified in 
Eq. (jni- 

B Lowest-order and one-loop amplitudes in terms of standard matrix ele- 
ments 

Both lowest-order and one-loop amplitudes can be written as linear combination of 
the SME defined in Eq. (JQ). The contribut ions to the lowest-order amplitude MIq, as 
defined in Eq. ()2.16|) . take the simple form 

^" " k (s - Mi)(s,, - Mir ' 



\AVt,<j ^ 3 9vc 9vt 9m 



X 



+ A^r + 2A^ij - MMr - Mr 

Ml'- = E 



,Vt,a _ „3 9 ve 9wt 9m 

t2\( ^ ^2\ 



r=± (s-M2)(s23-m2; 

X [Mr^ + Ml'-^ - 2Mg'-^ - mt(7Wr - Mri] ■ (B.l) 
The one-loop amplitude M.^ is much too lengthy to be written down; it is of the form 

28 

-Ml = E E F^^^M^, (B.2) 

i=l T=± 

where the coefficient functions Fr depend on couplings and scalar products of external 
momenta and contain all loop integrals. Note that here r is related to the chirality 
projectors u-r, whereas in Section |23 ti, r2 denote the helicities of the external fermions. 
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